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Organic-inorganic hybrid halide perovskites, in which the A cations of an ABX3 perovskite are re¬ 
placed by organic cations, may be used for photovoltaic and solar thermoelectric applications. In this 
contribution, we systematically study three lead-free hybrid perovskites, i.e., methylammonium tin 
iodide CHaNHsSnls, ammonium tin iodide NFLtSnU, and formamidnium tin iodide HC(NH2)2Snl3, 
by first-principles calculations. We find that in addition to the commonly known motif in which the 
corner-shared Snlg octahedra form a three-dimensional network, these materials may also favor a 
two-dimensional (layered) motif formed by alternating layers of the Snl6 octahedra and the organic 
cations. These two motifs are nearly equal in free energy and are separated by low barriers. These 
layered structures features many flat electronic bands near the band edges, making their electronic 
structures significantly different from those of the structural phases composed of three-dimension 
networks of Snl6 octahedra. Furthermore, because the electronic structures of HC(NH2)2Snl3 are 
found to be rather similar to those of CH3NH3S11I3, formamidnium tin iodide may also be promising 
for the applications of methylammonium tin iodide. 

PACS numbers: 61.50.-f, 31.15.E-, 72.20.Pa 


I. INTRODUCTION 

Ternary compounds with the ABX3 cubic perovskite 
structure, or simply perovskites, host an enormous num¬ 
ber of functionalities, e.g., ferroelectricity, colossal mag¬ 
netoresistance, and high thermopower. Interest in this 
family has further been fuelled by the development 
of organic-inorganic hybrid prerovskites in which A is 
replaced by an organic cation like methylammonium 
CH3NH3 or formamidinium HC(NH 2 ) 2 PKi3 While their 
photovoltaic applications with an efficiency of > 20% 11 
capture most of the recent attentionP®21 some hybrid 
halide perovskites, e.g., CH3NH3PM3 and CH3NH3S11I3, 
were also suggested to have large Seebeck coefficientP 0 
This implies that they may also be applicable for solar 
thermoelectric generators,® 0 converting the concent rated 
sunlight into electricity using the Seebeck effect 
These compelling applications promote a great deal of 
studies on this new class of materials, many of them are 
computation-based. 20 33 33 

Physical properties and hence, the efficiency of hybrid 
perovskites in such the applications are sensitive with 
their material structure, which turns out to be fairly 
complicated, as will be shown in this work. In the sim¬ 
ple cubic ABX 3 perovskite structure, A cations sit at 
the centers of the cages formed by a three-dimensional 
(3D) network of connected BX 6 octahedra. When large, 
anisotropic and polar organic cations are introduced at 
these sites, this structure is dramatically deformed, re¬ 
sulting two new major structural motifs. The first mo¬ 
tif, which is characterized by strongly distorted 3D net¬ 
work of BX 6 oct ahedra, has widely been examined by 

computationsPsEMS!] j n the sec0 nd motif, the 3D net¬ 
work is completely broken along one dimension, forming 
alternating layers of connected BXg octahedra and the 
organic cations. Consequently, structures of this two¬ 


dimensional (2D) layered motif are significantly different 
from the 3D structures. 

Layered st ructures in hybrid perovskites are 
ubiquitous ED3HH1 but computational effort devoted 
to them is surprisingly limited. In a work based on 
first-principles calculations,^ the experimentally re¬ 
solved structure of C4H9NH3PM4 was found to have 
a slightly wide band gap E g accompanied with flat 
electronic band dispersion nearby. These characteristics 
were also suggested^ for several nearly-2D structures 
of CH 3 NH3PbBr3 and CH 3 NH 3 PbCl3, signaled by the 
calculated long/short/short/long pattern of Pb-X axial 
distances (X = Br, Cl). Different from the 2D struc¬ 
ture examined for C4HgNH3Pbl4p^ the 2D structures 
predicted^ for CH 3 NH 3 PbBr3 and CH 3 NH 3 PbCl3 
share almost all the essential features of the 3D mo¬ 
tif, except slightly longer Pb-X axial distances. In a 
series of hybrid perovskites with rather large organic 
cations, e.g., [NH 2 C(I=NH 2 Y|2(CH3NH3) m Sn m I 3m+2 g 
[(CH 3 ) 3 NCH 2 CH 2 NH 3 ]Snl451(C 6 H5C 2 H4NH3) 2 Snl4pl 
and [CH 3 (CH 2 )iiNH3]SnI 3 f 17 the layered motif is unam¬ 
biguously dominant. Presumably, the 3D motif can not 
be preserved when such the large cations are introduced 
into the lattice. Even if the X anions are substituted by 
polyanions like BFj/ or PF4 , layered structures have also 
been suggested and studied computationally, revealing 
many flat electronic bands near the band edges P^ 

In this paper, we focus on three organic-inorganic 
hybrid lead-free halide perovskites, namely methylam¬ 
monium tin iodide CHsNFLSnL, ammonium tin iodide 
NELiSnLj, and formamidnium tin iodide HC(NH 2 ) 2 Snl3 
at the level of density functional theory (DFT)P^EH 
By employing the minima-hopping structure prediction 
method- 36 3 - we predict many low-energy structures for 
these hybrid perovskites, each of them exhibits ether the 
2D or the 3D motif. For each material, the lowest- 
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FIG. 1 . (Color online) Lowest-energy structures of 
CH 3 NH 3 SnI 3 , NH 4 SnI 3 , and HC(NH 2 ) 2 SnI 3 : the 3 D struc¬ 
ture on the left and the 2 D structure on the right. Sn, I, C, 
N, and H atoms are shown in yellow, purple, brown, lavender, 
and pink colors, respectively. 


energy 2D and 3D structures are separated by a low 
energy barrier, roughly 30 — 40 meV per atom. On 
the other hand, the 2D structures are characterized by 
flat electronic bands near the band edges, clearly dif¬ 
ferentiating themselves from the conventional 3D struc¬ 
tures. This observation suggests that the layered struc¬ 
tures of hybrid perovskites deserve further comprehensive 
investigations. We also find that the electronic struc¬ 
tures of HC(NH 2 ) 2 Snl 3 are rather similar to those of 
CHsNHsSnIa, suggesting formamidnium tin iodide may 
also be applicable for the photovoltaic and thermoelectric 
applicatons of methylammonium tin iodide. 


II. COMPUTATIONAL METHODS 

Our DFT calculations were performed with VASP 
package,®® employing the projector augmented wave 
(PAW) formalism.®] A basis set of plane waves with 
kinetic energies up to 500 eV was used to represent 
the wavefunctions while the local density approximation 
(LDA) functional was used for the exchange-correlation 
(XC) energies. The possible relativistic effects from the 
spin-orbit coupling (SOC) were included in our calcu¬ 
lations. The Brillouin zone of the examined structures 
was sampled by Monkhorst-Pacd® k-point meshes of 
7x7x7. The equilibration of the examined structures 
was assumed when the residual atomic forces exerting on 
the ions are smaller than 10 -2 eV/A. 

Low-energy structures of the hybrid perovskites were 
predicted using the minima-hopping method. 36 ]®] Typi¬ 
cally, a configurational space search is started by locally 


TABLE I. Calculated DFT energies Edft of the 2 D structures 
of three hybrid perovskites, given in meV/atom, with respect 
to those of the 3 D structures. Results were obtained by using 
various numerical prescriptions described in the text. 


Material 

LDA+SOC LDA PBE HSE 06 vdW-DF 2 

CH 3 NH 3 SnI 3 

7 

6 

10 

7 

2 

NH 4 SnI 3 

-6 

-7 

0 

-2 

-15 

HC(NH 2 ) 2 SnI 3 

1 

-1 

1 

-1 

-4 


optimizing an initial structure by a DFT code (VASP 
in this work). From the optimized structure, the en¬ 
ergy landscape is then constructed from the DFT en¬ 
ergy -Edft and explored by multiple attempts, each of 
them consists of a molecular-dynamics run for escaping 
the current local minimum and a local geometry opti¬ 
mization run on the obtained structure. The molecular- 
dynamic trajectory is chosen along the low-curvature di¬ 
rections of the energy landscape for improving the predic¬ 
tion efficiency, as suggested by Bell-Evans-Polanyi princi¬ 
ple. This is one of several feedback mechanisms and tech¬ 
niques that are employed in the minima-hopping method. 
This method has been successfully used for materials of 
various classes, including inorganic,®®! organic,®! and 
inorganic-organic hybrid compounds .®® 2 


III. LOW-ENERGY STRUCTURES 

Many low-energy structures of CH 3 NH 3 S 11 I 3 , 
NH 4 SnI 3 , and HC(NH 2 ) 2 SnI 3 were predicted, exhibiting 
either the 3D motif or the 2D motif. In this work, we 
focus on six lowest-energy structures, one 2D and one 3D 
structure for each of the three materials. In Fig. [TJ these 
six structures are visualized. Obviously, the topology of 
the ideal cubic perovskite structures is preserved in the 
structures of 3D motif, keeping them in a pseudo-cubic 
geometry, as experimentally described .^ 0 On the other 
hand, layers of BXg octahedra are unambiguously 
formed in the 2D structures and strongly shifted with 
respect to each other. Moreover, the simulated x-ray 
diffraction (XRD) shown in Fig. 1 of the Supplemental 
Material,® indicates that the examined 3D structures 
of CH 3 NH 3 Snl 3 and HC(NH 2 ) 2 Snl 3 are essentially 
identical with the pseudo-cubic a phases experimentally 
resolved® for these perovskites and computationally 
studied elsewhere. 2 ® In Sec. [vj this conclusion will 
further be confirmed by the excellent consistence be¬ 
tween the calculated and the measured band gaps of 
CH 3 NH 3 Snl 3 and HC(NH 2 ) 2 Snl 3 . The crystallographic 
information and simulated XRD patterns of all the 2D 
and 3D structures can also be found in the Supplemental 
Material.® The crystallographic information is also 
reported in http://khazana.uconn.edu. 

For any hybrid perovskite, the lowest-energy 2D and 
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FIG. 2. (Color online) Low-frequency regions of the cal¬ 
culated phonon band structures of the thermodynamically 
most stable 2D and 3D structures predicted for CHsNHsSnD, 
NH 4 S 11 I 3 , and HC(NH 2 ) 2 Snl 3 (the full vibrational frequency 
spectra of these structures span up to ~ 100 THz). For conve¬ 
nience, phonon bands with imaginary frequencies are shown 
as those with negative frequencies. Red dashed lines sketch 
the numerical error of ~ 0.3 THz typically encountered by 
phonon calculations. 


3D structures are energetically competing with rather 
small energy difference. In case of NH 4 S 11 I 3 , the 2D struc¬ 
ture is lower than the 3D structure by ~ 6 meV/atom. 
On the other hand, the 2D structure of CHaNHsSnlg 
and HC(NH 2 ) 2 Snl 3 is higher than the 3D structures 
by roughly 7 meV/atom and 1 meV/atom, respectively. 
To further confirm the thermodynamic stability of the 
2D structures with respect to the 3D structures, we 
performed additional calculations for the DFT energy 
Soft using different numerical schemes, employing the 
LDA, the Perdew-Burke-Ernzerhof (P BE)p^ and the 
hybrid Heyd-Scuseria-Ernzerhof (HSE06) 55 * 56 ! XC func¬ 
tionals (without SOC). To somehow capture the long- 
range dispersion interaction, we also used a prescrip¬ 
tion which incorporates the vdW-DF2 non-local density 
functional!®^ The results are shown in Table [IJ indicating 
that the 2D and 3D motifs of HC(NH 2 ) 2 Snl 3 are essen¬ 
tially the same in energy. For NH 4 Snl 3 , the 2D struc¬ 
ture is slightly more stable than the 3D structure while 
for CHsNHsSnIs, the 3D structure is slightly more stable 



FIG. 3. (Color online) Helmholtz free energy of the energeti¬ 
cally most favorable 2D structures of CHsNHsSnD, NHiSnD, 
and HC(NH 2 ) 2 Snl 3 , given with respect to that of the most 
stable 3D structure of the same material (dotted line). 


than its 2D counterpart. For all three hybrid perovskites, 
the 2D and 3D structures are essentially comparable in 
energy. 

The phonon frequency spectrum of a theoretically pre¬ 
dicted structure is typically used for examining its dy¬ 
namical stability. We used the su perce ll approach, as 
implemented in PHONOPY package, ^ — ^ to perform the 
relevant frozen-phonon calculations for all the six struc¬ 
tures structures examined. For each of them, a 2 x 2 x 2 
supercell was used for the phonon calculations. We note 
that calculations of this type are performed within the 
harmonic approximation on the structures which were 
determined at zero temperature. The low-frequency re¬ 
gion (5 THz and below) of the calculated phonon spectra 
is shown in Fig. [ 2 ] while the full spectra (spanning up to 
~ 100 THz) are given in the Supplemental Material . 53 
Fig. H indicates that within the known numerical er¬ 
ror of phonon calculations (~ 0.3 THz, or roughly 0.1 
meV/atom for our structures) j 47 l 6 ° l the six structures in¬ 
deed correspond to the local minima of the energy land¬ 
scape, and hence, are dynamically stable. 

At finite temperatures T, the thermodynamic stabil¬ 
ity of these structures may be accessed by estimating the 
Helmholtz free energy F(T) = £45 FT + £ vib (T).®® In 
this expression, A v jb(T) is the entropic contribution to 
F(T) from the lattice vibrations. This term can be esti¬ 
mated by the lattice dynamics approach within the har¬ 
monic approximation from the calculated phonon density 
of states g(nj) atPSHIl 


F v ib{T) = 3NknT / dujg(oj) In 


2 sinh 



(!) 

Here, fce is the Boltzmann constant, N the number of 
degrees of freedom, and u> the phonon frequency. The 
phonon density of state g(ui) used for Eq. [I] were evalu¬ 
ated by sampling the Brilluion zone by a q-point mesh of 
9x9x9. We show in Fig. [ 3 ] the free energy F(T) calcu- 
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FIG. 4 . (Color online) Computed minimum energy path¬ 
ways from the 3 D (left) to the 2 D (right) structures of 
CH3NH3S11I3, NH 4 SnI 3 , and HC(NH 2 )2SnI 3 . Curves were 
interpolated from the calculated energies of the images, given 
by symbols. 


lated for the six structures examined, indicating that the 
2D structure of NF^Snlg is thermodynamically more sta¬ 
ble than the 3D structure at low temperatures but grad¬ 
ually become less stable at T ~ 350K and above. In cases 
of CH 3 NH 3 Snl 3 and HC(NH 2 ) 2 SnI 3 , the 3D structure is 
always favorable over the 2D structure within the tem¬ 
perature range examined. Overall, in all three materials, 
the energy differences between the 2D and 3D structures 
are relatively small and comparable with k^T. which is 
about 25 meV/atom at room temperature. 


IV. ACTIVATION BARRIERS 




rxMr RMXRrxMr rmxr 


FIG. 5 . (Color online) Computed electronic bands of 
the 3 D and 2 D structures of CH3NH3S11I3, NFUSnU, and 
HC(NH2)2Snl3. Fermi energies are set to zero. 


cations, and then forming new Sn-I bonds to construct 
the 2D motif. 0 ' 3 ! Moreover, all of the organic cations are 
rigid and highly anisostropic, so their rotations will not 
be smooth. Consequently, the transformation from a 3D 
to a 2D structure generally encounters a sequence of bar¬ 
riers, as can be seen in Fig. [4] At operating temperatures 
(room temperature and above), these barriers are how¬ 
ever small, being comparable with k^T. 


Because the 2D and 3D motifs of CH 3 NH 3 SnI 3 , 
NH 4 SnI 3 , and HC(NH 2 ) 2 SnI 3 are nearly equal in energy, 
it is interesting to examine the energy barriers between 
them. To find this information, we employed the solid- 
state (climbing image) nudged elastic band methocP^MI 
for estimating the minimum energy pathways between 
the 3D and the 2D structures of the hybrid perovskites. 
The calculated results are shown in Fig. [4j demonstrat¬ 
ing that for these materials, the 3D and 2D motifs are 
separated by low energy barriers, falling between 30 — 40 
meV/atom. The structural transition from a 3D to a 2D 
structure involves a series of steps, including breaking the 
long Sn-I “bonds” along the out-of-plane direction, rotat¬ 
ing both the in-plane Sn-I bonds and the whole organic 


V. ELECTRONIC STRUCTURES 

We performed electronic structure calculations for the 
2D and 3D structures of the examined hybrid perovskites, 
and the results are shown in Fig. [5j Interestingly, while 
the 3D structures are characterized by high-curvature 
parabolic bands, the 2D structures feature many flat 
bands next to the conduction band minimum (CBM) 
and the v alenc e band maximum (VBM) as previously 
addressed, 18 E31thereby leading to some high-density elec¬ 
tronic levels near the band edges. At the LDA+SOC level 
of DFT computations, the band gap E s of the 2D struc¬ 
tures is found (see Table 0 to be considerably higher 
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TABLE II. Calculated band gap of the 2D and 3D structures 
of three hybrid perovskites, given in eV. Results were obtained 
by using the HSE06 functional for the exchange-correlation 
energies while those obtained with LDA+SOC are given in 
parentheses for completeness. 


Material 2D 3D Literature 

CH 3 NH 3 SnI 3 3.00 (1.3) 1.22 (0.2) 1.20 - 1.3^ 
NH 4 SnI 3 3.05 (1.9) 1.65 (0.4) 
HC(NH 2 ) 2 SnI 3 2.75 (2.0) 1.20 (0.2) 1.4P1 


than that of the 3D counterparts while the CBM degen¬ 
eracy is unambiguously lifted. 

As the intrinsic band gap underestimation of DFT 
with (semi)local XC functionals like LDA and PBE is 
inevitable,^ we performed additional calculations for the 
electronic structures with the hybrid HSE06 XC func¬ 
tional. The obtained electronic band structures are given 
in the Supplemental Material while the calculated band 
gaps are shown in Table [TlJ For the 3D structure of 
CH 3 NH 3 SnI 3 and HC(NH 2 ) 2 SnI 3 , the calculated band 
gaps are 1.22 eV and 1.20 eV. These values are well con¬ 
sistent with the experimentally determined band gaps 
of 1.20 eV—1.35 eV and 1.41 eV, respectively ? 0 At the 
HSE06 level of DFT, the calculated band gap of the 2D 
structures is also significantly larger that that of the 3D 
structures, as revealed by calculations at the LDA+SOC 
level. Moreover, the electronic structures of the 3D struc¬ 
tures of CH 3 NH 3 SnI 3 and HC(NH 2 ) 2 SnI 3 (those exper¬ 
imentally realized) are rather similar, suggesting that 
HC(NH 2 ) 2 SnI 3 may also be promising for the photo¬ 
voltaic and thermoelectric applications of CH 3 NH 3 SnI 3 . 


VI. SUMMARY 

In summary, we have systematically studied three 
organic-inorganic hybrid halide perovskites, i.e., 
CH 3 NH 3 SnI 3 , NELjSnhj, and HC(NH 2 ) 2 SnI 3 , by first- 
principles calculations. We find that in addition to the 
3D network of BX 6 octahedra, the 2D motif formed by 
alternating layers of BXg octahedra and organic cations 
may also be favorable. When the ionic cations are 
large, the layered structures would become dominant, 
as unambiguously observed in experiments?^® While 
these two structural motifs are energetically competing 
in CH 3 NH 3 SnI 3 , NILiSnLj, and HC(NH 2 ) 2 SnI 3 , they are 
separated by low barriers (~ 30 — 40 meV/atom). The 
2D layered structures are significantly different from the 
conventional 3D structures in terms of the electronic 
structures. On the other hand, the electronic structures 
of the 3D motifs of CH 3 NH 3 SnI 3 and HC(NH 2 ) 2 SnI 3 
are quite similar. In regard of the photovoltaic and ther¬ 
moelectric applications of the hybrid halide perovskites, 
the present results suggest that further in-depth inves¬ 
tigations of formamidnium tin iodide HC(NH 2 ) 2 SnI 3 
and the layered structures of the hybrid perovskites are 
desirable. 
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